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Advances in cryogenic sciences have enabled several observations of new low-temperature phys¬ 
ical phenomena including superconductivity, superfluidity, and Bose-Einstein condensates. Heat 
transfer is also critical in numerous applications including thermal management within integrated 
microelectronics and the regulation of plant-growth and development. Here we demonstrate that 
single-beam laser-trapping can be used to induce and quantify the local refrigeration of aqueous 
media through analysis of the cold Brownian dynamics of individual Yb®^-doped yttrium lithium 
fluoride (YLF) crystals in an inhomogeneous temperature field via forward light scattering and 
back-focal-plane interferometry. A tunable, NIR continuous-wave laser is used to optically trap 
individual YLF crystals with an irradiance on the order of 1 MW/cm^. Heat is transported out of 
the crystal lattice (across the solid / liquid interface) by anti-Stokes photoluminescence following 
upconversion of Yb®"*" excited states mediated by optical-phonon absorption. The cold Brownian 
motion (CBM) analysis of individual YLF crystals indicates local cooling by >21°C below ambient 
conditions suggesting a range of potential future applications. 


I. INTRODUCTION 

In 1929, Pringsheim[l] proposed that solid-state ma¬ 
terials could experience refrigeration if they exhib¬ 
ited biased emission of anti-Stokes (blue-shifted) radi¬ 
ation relative to a fixed optical excitation wavelength. 
Epstein[2] and colleagues experimentally demonstrated 
this concept first in 1995 using rare-earth-doped fluo¬ 
ride glass materials (ZBLAN). More recently, it has been 
shown that solid-state rare-earth-doped yttrium-lithium- 
fluoride (Yb^+:YLF) crystals can be cooled to cryogenic 
temperatures [3] in vacuo using a continuous-wave NIR 
laser excitation. Furthermore, the laser refrigeration 
of doped yttrium-aluminum-garnet (Yb^+iYAG) mate¬ 
rials has recently been reported in air at atmospheric 
pressure [4]. Anti-Stokes photoluminescence has also been 
shown[5] to cool cadmium sulfide (CdS) nanoribbons in 
vacuo by as much as 40 °C below room temperature. 

In contrast to anti-Stokes processes, optomechanical 
laser refrigeration has also been demonstrated based on a 
novel mechanism of angular momentum transfer between 
a circularly polarized laser and a birefringent crystal[6]. 
To date, laser-refrigeration of nanocrystals in aqueous 


* These authors contributed equally to this work. 
"I" Corresponding author: peterpz@uw.edu 


media has not been reported stemming primarily from 
the large near-infrared (NIR) optical absorption coef¬ 
ficient of water [7]. It has remained an open question 
whether these known cooling materials could act to re¬ 
frigerate aqueous media and undergo hypothesized cold 
Brownian motion (CBM) [8, 9]. In this work, we demon¬ 
strate the local refrigeration of aqueous media by > 21 °C 
below ambient conditions following anti-Stokes photolu¬ 
minescence from optically trapped[10], rare-earth-doped 
YLF nanocrystals undergoing CBM. 


II. RESULTS & DISCUSSION 
A. Nanocrystal Synthesis & Characterization 

Pioneering efforts to cool Yb^+:YLF materials in vacuo 
have relied on the growth of high-purity YLF single¬ 
crystals using an air- and moisture- free Czochralski 
process [11]. In the experiments reported here, a low- 
cost modified hydrothermal synthesis[12] of Yb^+iYLF is 
used to prepare YLF crystals shown in Fig. 1. Scanning 
electron microscopy reveals that YLF crystals exhibit 
a truncated tetragonal bipyramidal (TTB) morphology 
(Fig. lb). X-ray diffraction shows that the YLF crys¬ 
tal has a Scheelite structure (Fig. Ic). Bright field / 
HAADF TEM imaging (Fig. Id/e) and electron diffrac- 
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FIG. 1. Synthesis and characterization of YLF crystals, a, Schematic of Scheelite crystal structure of YLF with In/a 
space group symmetry, b, Scanning electron microscope image of a faceted (Yb®^)o.i(Y^+)o. 9 LiF 4 particle exhibiting TTB 
morphology. Scale bar = 1 /rm. c, Powder x-ray diffraction pattern of YLF crystals following hydrothermal synthesis indicating 
a pure Scheelite crystal phase. Inset: schematic of TTB morphology relative to YLF’s unit cell, d, Bright field transmission 
electron microscope (TEM) image of an individual Yb^^iYLF grain; scale bar = 200 nm. Inset: high-resolution TEM image 
taken from the indicated region; scale bar = 2 nm. e, High-angle annular-dark-field (HAADF) image of the YLF grain in panel 
b showing regions of high contrast suggesting the presence of polycrystalline domains. Inset: select area electron diffraction 
from the indicated region, f, X-ray fluorescence compositional-analysis-spectrum of an individual YLF crystal taken within the 
TEM confirming the elemental crystalline composition including Y, Yb, and F species. 


tion suggest that the TTB materials are polycrystalline 
and likely form through an oriented attachment [13] pro¬ 
cess of nanocrystalline grains (Fig. le, inset). 


B. Optical Trap Setup 

A laser trapping system shown in Fig. 2 was 
used to optically trap individual Yb^+:YLF crystals. 
Briefly, a single-beam laser trap was used to extract 
the surrounding local temperature profile of optically 
trapped rare-earth-doped YLF particles by analysis of 
calibrated power spectra of the forward-scattered laser 
radiation[14], which is used to quantify the CBM of the 
particle. In order to minimize fluid absorption and heat¬ 
ing at the trapping wavelengths[15], laser trapping ex¬ 
periments were performed in D 2 O unless explicitly stated 
otherwise. The laser trap setup is outlined in Materials 
& Methods, and the CBM temperature analysis is de¬ 
scribed in Appendix VIA. 


C. CBM of Locally Refrigerated Nanocrystals 

A bright-held micrograph for a characteristic optically- 
trapped Yb^“'':YLF crystal is shown in Fig. 3a. The 
dependence of laser refrigeration on the trapping lasers 


pump wavelength is shown in Fig. 3c, where YLF crystals 
doped with 10% Yb3+ are observed to cool from 19 °C 
at a 5.9 MW/cm^ trapping irradiance to 4 °C at a 25.5 
MW/cm^ trapping irradiance when trapped with A = 
1020 nm, which is resonant with ytterbium’s E4-E5 tran¬ 
sition shown in Fig. 3b. The same Yb^“'':YLF crystals 
are shown to heat from 40 °C to 47 °C when trapped at 
the same respective irradiances with A = 1064 nm, which 
is energetically insufficient to pump the E4-E5 resonance 
and subsequently cannot initiate upconversion-mediated 
cooling. As a control, 1.01 /rm diameter amorphous sil¬ 
ica microspheres are observed to heat from 39 °C to 46 
°C when trapped with the A = 1020 nm trapping wave¬ 
length. 


D. Ratiometric Thermometry of Codoped 
Nanocrystals 

Furthermore, it has been recently reported that strong 
visible upconversion in rare-earth codoped nanocrys¬ 
tals can be used for efficient biological imaging and 
labeling[16]. The long (ms) lifetime of photoluminescence 
from rare-earth point defects may also be used to infer 
temperature changes through ratiometric thermometry 
by analysis of Boltzmann thermal populations[17], which 
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FIG. 2. Schematic of laser trapping instrument. An optically trapped YLF crystal in an aqueous fluid chamber. A 
piezostage driven at 32 Hz produces a peak in the quadrant photodiode (QPD) power spectrum which is used to extract 
a calibrated diffusion constant. The particle’s temperature (Tp) and local temperature profile is then extracted using cold 
Brownian motion analysis. 


is given by the equation: 

h 

— oc exp 

h 


/ -(E2-E^) \ 
\ hT ) 
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In brief, changes in the ratio of the integrated emis¬ 
sion bands I 2 and Ii that stem from transitions be¬ 
tween energy states E 2 and Ei, respectively, and a com¬ 
mon ground state are directly correlated to a change 
in the particle’s temperature. Codoping YLF crystals 
with both Yb^+ and Er^+ ions leads to a thermalized 
Boltzmann distribution between the E 2 (^Hii/ 2 ) and Ei 
(^ 83 / 2 ) manifolds of Er^"*" and an intense green upconver- 
sion emission that is visible to the unaided eye, as shown 
in Fig. 4a. 

Utilizing photoluminescence spectroscopy with opti¬ 
cally trapped codoped upconverting nanocrystals pro¬ 
vides a unique capability of testing the particle-to- 
particle variability in emission. For the codoped 2%Er^+, 
10%Yb^+:YLF particles reported, substantial fluctua¬ 
tions in upconversion photoluminescence was observed, 
making bulk calibration measurements inadequate for in¬ 
dividual ratiometric temperature measurements as de¬ 
tailed in Appendix VIC. Thus, ratiometric thermome¬ 
try can be used only qualitatively during laser trapping 
experiments to observe changes in temperature as the 
trapping irradiance is increased, as shown in Fig. 4b/c. 
The increase (decrease) in the logarithmic ratio of I 2 to 
Ii (Fig. 4b/c) with increasing irradiance reflects an in¬ 
crease (decrease) in the internal lattice temperature[18]. 


which agrees well with the observed temperature changes 
measured via laser trapping light scattering temperature 
analysis (Fig. 4d). Specifically, laser trapping analysis 
of the particles’ CBM indicates that codoped 2%Er^+, 
10%Yb^''':YLF undergoes laser refrigeration (AT = -4.9 
± 2.8 °C) when trapped at A = 1020 nm and heating 
(AT = 21.8 ± 10.11 °C) when trapped at A = 975 
nm. Furthermore, it has been proposed recently that 
codoping YLF crystals with other upconverting rare- 
earth ions can enhance cooling through energy transfer 
enhanced cooling[19]. The laser trapping and refrigera¬ 
tion of 10%Yb^“*':YLF nanocrystals at A = 1020 nm in bi¬ 
ological buffers such as phosphate-buffered saline (PBS) 
and Dulbecco’s Modified Eagle Medium (DMEM) is also 
possible, as demonstrated in Table 1. 


TABLE I. Local cooling of Yb^"'"-doped YLF crystals 
in various media. 


Solvent 

AT = (Tp - To) (°C) 

AT Std. Dev. (°C) 

D 2 O 

-15.0 

4.1 

DI Water 

-14.7 

3.8 

PBS 

-14.9 

4.3 

DMEM 

-11.2 

6.3 


. * Temperature changes are determined between the laser irradi- 
ances 5.9 and 25.5 MW/cm^, where To = 25 '^C. 
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FIG. 3. Laser refrigeration of optically trapped YLF 
microcrystals. a, Optical micrograph of an optically 
trapped YLF crystal; scale bar = 3 /rm. b, Crystal field 
energy level confignration of Yb^^ dopant ions and em¬ 
ployed cooling scheme, c, Extracted temperature of opti¬ 
cally trapped particles in D 2 O. Yb^+-doped YLF particles 
are shown to cool when trapping wavelength is resonant with 
the E4-E5 transition {Xtrap = 1020 nm) but heat when the 
trapping wavelength is below the transition {Xtrap ~ 1064 
nm). As a control, amorphous silica microspheres heat when 
trapped at Xtrap = 1020 nm. 


E. Outlook 

These results illustrate the potential of using singly- 
and co-doped YLF nanocrystals as a platform for physi¬ 
ological refrigeration, biological imaging, and in situ ra- 
tiometric thermometry. Potential applications for these 
materials include triggering and probing fundamental 
temperature-controlled biological processes[20, 21]. The 
ability to measure and to control temperature could en¬ 
able this technology to investigate kinetics and temper¬ 
ature sensitivity of basic cellular processes, including 
ion channels [ 22 ], or to act as simultaneous theranostic- 
hypothermia agents to identify and treat cancerous 
tissues [23]. 

Analyzing the CBM of a nanocrystal dispersed in a 
liquid phase to measure the nanocrystal’s temperature 
also provides the unique capability to extract the lo¬ 
cal temperature gradient in the medium surrounding the 
trapped nanocrystal. Since the aspect ratio of the TTB 
materials used here is near unity, we approximate the ra¬ 
dius R of the particles using an equivalent-sphere model 
and can extract the local temperature field a distance 
r from the particles’ surface (at temperature Tp, exclud¬ 
ing the temperature discontinuity at the particle’s surface 


from the Kapitza resistance[24]), which is given by[25]: 

T{r)=To + ^{Tp-To) (2) 

where Tq is the bath temperature of the medium. 
Given that the average radius of the Yb^+:YLF parti¬ 
cles trapped at A = 1020 nm in Fig. 3c is Ra«g = 764 
± 293 nm, Tq = 25 °C, and Tp^avg = 3.4 °C at 25.5 
MW/cm^ irradiance, the distance away from the par¬ 
ticle where the temperature increases to within 1 % of 
Tq is 6.9 /rm (Fig. 2). Future synthetic developments 
with YLF host crystals could enhance the resonant op¬ 
tical absorption of Yb^"*" through morphology dependent 
resonances[26], or reduce unwanted background impurity 
absorption to a level capable of locally freezing water at 
atmospheric or elevated pressures. 


III. MATERIALS & METHODS 

A. YLE Synthesis 

The following synthesis was performed following mod¬ 
ifications to Lu, C. et. al[12]. Yttrium oxide (Y 2 O 3 ), 
ytterbium oxide (Yb 203 ) and erbium oxide (Er 2 03 ) are 
of 99.99% purity and used as purchased from Sigma- 
Aldrich. Yttrium nitrate (Y(N 03 ) 3 ), ytterbium nitrate 
(Yb(N 03 ) 3 ) and erbium nitrate (Er(N 03 ) 3 ) are obtained 
by dissolving the oxide in concentrated nitric acid at 60 
°C while stirring for several hours until excess nitric acid 
is removed. The residual solid is then dissolved in Milli- 
pore deionized (DI) water to achieve a stock concentra¬ 
tion of the respective nitrate. Lithium fluoride (LiF), ni¬ 
tric acid (HNO 3 ), ammonium bifluoride (NH 4 HF 2 ) and 
ethylenediaminetetraacetic acid (EDTA) are analytical 
grade and used directly in the synthesis without any pu¬ 
rification. The following preparation uses the synthesis of 
2 %Er^+, 10 %Yb^“'':YLiF 4 as an example. 7.04ml of 0.5M 
Y(N 03 ) 3 , 0.8ml of 0.5M Yb(N 03)3 and 0.16ml of 0.5M 
Er(N 03)3 are mixed with 1.17g EDTA in 5ml Millipore 
DI water at 80 °C while stirring for Ih. This is solution 
A. Subsequently, 0.21g of LiF and 0.68g of NFI 4 HF 2 are 
dissolved in 7ml Millipore DI water at 70 °C while stir¬ 
ring for Ih to form solution B. Solutions A and B are 
mixed together while stirring for 20 min to form a ho¬ 
mogeneous white suspension which is then transferred to 
a 23ml Teflon-lined autoclave and heated to 220 °C for 
72h. After the autoclave cools to room temperature, the 
2 %Er^+, 10 %Yb^''":YLiF 4 particles can be recovered by 
centrifuging and washing with ethanol and Millipore DI 
water three times. The final white powder is obtained by 
calcining at 300 °C for 2h. 10 %Yb^+:YLiF 4 particles are 
achieved using the same method. 
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FIG. 4. Upconversion and ratiometric thermometry of codoped YLF. a, Bright-field optical micrograph showing a 
codoped 2%Er^"'",10%Yb®"'":YLF particle in Brownian motion (top-left) and a dark-field optical micrograph of the crystal when 
trapped with Xtrap = 1020 nm (bottom-left). Scalebar = 4/im. Upconverted photoluminescence can be seen with the unaided 
eye (right), b, Photoluminescence spectra of the corresponding dark-field image showing the integration regions I 2 and Ii, 
representing emission from Er^+ energy states E 2 (^Hii/ 2 ) and Ei ('‘S 3 / 2 ) to the ground state Eground (‘* 115 / 2 ), respectively, 
c, Natural logarithm of the ratio I 2 /I 1 showing a linear increase (top) with laser irradiance at Xtrap = 975 nm and a linear 
decrease (bottom) with laser irradiance at Xtrap = 1020 nm. d, Laser refrigeration of the codoped YLF crystal analyzed in c 
measured via cold Brownian motion analysis. 


B. TEM Characterization 

Bright field and STEM HAADF images were taken on 
a FEI Tecnai G2 F20 at an accelerating voltage of 200 
keV. Select area electron diffraction (SAED) images were 
taken with a camera length of 490 mm. EDS spectra were 
obtained with a 60 second acquisition time. The spectra 
were then processed by subtracting the background and 
smoothing the peaks. 


C. SEM Characterization 

Secondary electron images were taken on an FEI Sirion 
at an accelerating voltage of 5 keV. 


D. XRD Characterization 

Powder x-ray diffraction (XRD) patterns are obtained 
on a Bruker F8 Focus Powder XRD with Cu K (40kV, 
40mA) irradiation (A=0.154nm). The 26 angle of the 
XRD spectra is from 10° to 70° and the scanning rate 
is 0.01°s“*. The one minor unlabeled peak in the XRD 
spectra at 26 = 44.9° is attributed to a small amount of 
unreacted LiF precursor ((200) peak). 


E. Laser Trapping Description 

The laser tweezer setup is a modihed modular opti¬ 
cal tweezer kit (Thorlabs, OTKB), where the original 
condenser lens has been replaced with a lOx Mitutoyo 
condenser (Plan Apo infinity-corrected long WD objec¬ 
tive, Stock No. 46-144). The lOOx objective focusing lens 
has a numerical aperture of 1.25 and a focal spot of 1.1 
/im. The quadrant photodiode and piezostage were inter¬ 
faced to the computer through a DAQ card (PCIe-6361 
X Series, National Instruments) and controlled through 
modified MATLAB software (Thorlabs). Experimental 
chambers were prepared as follows. Several microliters of 
the nanocrystal/aqueous medium dispersion were trans¬ 
ferred by a pipette into a chamber consisting of a glass 
slide and glass coverslip. The edges of the glass slide 
and the glass coverslip were then sealed with a 150-//m- 
thick adhesive spacer (SecureSeal Imaging Spacer, Grace 
Bio-labs). Nanocrystals were trapped at the center ( 75 
/im from the surface) of the temperature controlled per¬ 
fusion chamber (RC-31, Warner Instruments) and held 
at To = 25 °G while voltage traces were recorded at the 
quadrant photodiode (QPD) for 3 seconds at a sample 
rate of 100 kHz. The QPD voltage signal was calibrated 
by oscillating the piezostage at 32Hz and an amplitude 
of 150nm peak-to-peak during signal acquisition, as out¬ 
lined in reference [27]. Trapping data was acquired using a 
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diode-pumped solid state YbiYAG thin-disk tunable laser 
(VersaDisk 1030-10, Sahajanand Laser Technologies) at 
a wavelength of 1020 nm, a 975 nm pigtailed Fiber 
Bragg Grating (FBG) stabilized single-mode laser diode 
(PL980P330J, Thorlabs), as well as a solid-state Nd:YAG 
1064 nm (BL-106G, Spectra-Physics) at an irradiance of 
5.9, 10.7, 14.6, 21.2, and 25.5 MW/cm^. Each YLF cool¬ 
ing data point in Fig. 3c in the manuscript represents an 
average of 6 individual particles with an average radius 
of 764 nm with a standard deviation of 293 nm. Silica 
beads (SS04N/9857, Bangs Laboratories) were used for 
their monodisperse size distribution (1010 nm diameter), 
and they have shown to minimally heat when trapped 
with a laser tweezer at NIR wavelengths [28]. Electro¬ 
magnetic simulations of the interaction of the trapping 
laser with a YLF TTB were also performed to predict 
the stable trapping configurations of optically trapped 
YLF particles, detailed in Appendix VIE. Lastly, visible 
emission of Er^+ from Er/Yb codoped trapped YLF host 
crystals was detected using an Acton SpectraPro 500i 
spectrograph with a Princeton liquid-nitrogen cooled Si 
detector. 
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VI. APPENDICES 


A. Appendix I: Cold Brownian Motion 
Temperature Extraction with Laser Trapping 


Power spectra from the quadrant photodiode (QPD) 
voltage traces were processed according to Berg-Sprensen 
and Flyvbjerg [29] and used to calibrate the QPD traces 
following the method of Tolic-Nprrelykke et al. [27]. An 
experimental diffusion coefficient was then extracted by 
fitting the characteristic function for the experimental 


power spectra derived in Berg-Sprensen and Flyvbjerg 
[29]. Given that the temperature of the trapped particle 
is significantly different than the temperature sufficiently 
far from the laser focus, the particle-trap system is not 
isothermal and behaves according to nonequilibrium dy¬ 
namics. Thus, equating the experimental diffusion coef¬ 
ficient to nonisothermal Brownian dynamics necessitates 
the application of cold Brownian motion (CBM), as de¬ 
rived by Chakraborty et al. [8]. The CBM diffusion co¬ 
efficient is then related to the CBM temperature by: 


Dcbm 


kbTcBM 

1cbm{T) ’ 


( 3 ) 


where Dcbm is the CBM diffusion coefficient, kb is 
Boltzmann’s constant, Tqbm is the CBM temperature, 
and joBMiT) is the CBM Stokes drag. To the leading 
order of the temperature increment or decrement AT 
= (Tp - To), the temperature-dependence of the vis¬ 
cosity on Tcbm can be neglected, giving the effective 
temperature [8] 


Tcbm =To + j^^T. 


( 4 ) 


For a temperature dependence of the solvent viscosity 


of the form r](T) = r/oo exp 
drag is given by 


(T-Tvf) 


, the CBM Stokes 


7cbm(T) = QnRrjcBMiT), (5) 

where R is the particle radius, and rjcBMiT) is the 
temperature-dependent CBM viscosity that is related to 
the viscosity of the solvent at room temperature, tjq, by 


Vo 

Vcbm{T) 


I-b 


193 

4^ 



AT 

{To - Tvf) 




12563 2 

118098 



AT 

{To - Tvf) 


( 6 ) 


Equations 4-6 are then used in Eq. 3 to obtain Dcbm, 
which is subsequently compared to the experimental dif¬ 
fusion coefficient to determine the particle temperature 
Tp (excluding the temperature discontinuity at the par¬ 
ticle’s surface from the Kapitza resistance [24]). 

An alternative CBM temperature analysis using a 
semi-phenomenological expression for Dcbm that ap¬ 
proximately accounts for higher order terms in AT (Eq. 
15 of the supporting online materials of Chakraborty 
et al. [8]) yields consistent results, indicating that these 
higher order corrections are negligible, for our purposes. 
For the experiments reported here, the VFT viscosity pa¬ 
rameters were fit to experimental data and are as follows: 

D 2 O 

Voo = 3.456 * 10“® Pa-s 
A = 478.6 K 
Tvf = 160 K, 


and 


DI water, PBS, DMEM 

rjoo = 2.664 * 10“® Pa-s 
A = 536.5 K 
Tvf = 145.5 K. 


VFT viscosity parameters for DI water, PBS (O.OIM, 
pH 7.4; Sigma P5368), and DMEM (IX, high glucose, 
pyruvate; Life Technologies Cat. #11995-065) were as¬ 
sumed to be equivalent since it has been reported that 
water viscosity can be used for purposes of modeling par¬ 
ticle transport in non-serum containing media[30]. 



a 



FIG. 5. Laser trap temperature extraction error anal¬ 
ysis. A 1.01 fim diameter silica sphere trapped in D 2 O us¬ 
ing a trapping wavelength of 975 nm in a temperature con¬ 
trolled chamber held at To = 25 °C. QPD data was acquired 
10 consecutive times at an irradiance of 5.9 MW/cm^ and 
the temperature of the sphere was determined using (red 
circles, Tcbm) the cold/hot Brownian temperature extrac¬ 
tion method detailed in Appendix VIA, and (blue squares, 
Tciasaic) the classical isothermal assumption method. The 
error bars represent the standard deviation. 
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B. Appendix II: Laser Trap Temperature 
Extraction Error Analysis 


In order to determine the accuracy of the temper¬ 
ature extraction method presented here, silica beads 
(SS04N/9857, Bangs Laboratories) are used for their 
monodisperse size distribution (1010 nm diameter) and 
their known capacity to heat minimally when trapped 
with a NIR laser trap[28]. The same experiment that was 
performed on YLF particles in this study was performed 
on the silica beads in D 2 O using a trapping wavelength 
of 975 nm and an irradiance of 5.9 MW/cm^. For this 
experiment, a temperature controlled chamber was used 
(RC-31, Warner Instruments), with the chamber held at 
To = 25 °C while the same bead was trapped and QPD 
data acquired 10 consecutive times. The cold Brownian 
analysis outlined above in Appendix VIA was then ap¬ 
plied to the resulting QPD signals, and the temperature 
results are shown in red circles in Fig. 5. The average cal¬ 
culated particle temperature (red dashed line) is Tcbm 
= 35.9 °C with a standard deviation (red shaded area) 
of ± 6.8 °C. For comparison, the same data was also an¬ 
alyzed with the method assuming isothermal conditions, 
and the temperature results are shown in blue squares 
in Fig. 5. The average calculated particle temperature 
(blue dashed line) is Tclassic = 29.5 °C with a standard 
deviation (blue shaded area) of ±2.8 °C. 


FIG. 6. Ratiometric thermometry variability, a, Select 
photoluminescence spectra from the measurements in (b) of 
multiple optically trapped, codoped 2%Er®''',10%Yb^+:YLF 
particles showing the variability in the integration regions 
I 2 and Ii, representing emission from Er^^ energy states E 2 
(^Ilii/ 2 ) and El ("‘S 3 / 2 ) to the ground state Eground ('‘I 15 / 2 ), 
respectively, b, Natural logarithm of the ratio I 2 /I 1 of dif¬ 
ferent trapped particles showing the particle-to-particle vari¬ 
ation. 


C. Appendix III: Ratiometric Thermometry 
Calibration 


Calibration of the ratiometric thermometry on an in¬ 
dividual particle is non-trivial in that there exist fluc¬ 
tuations in particle morphology and dopant concentra¬ 
tion/distribution that affect the kinetics of the upcon- 
version efficiency and overall spectral profile of the Ei 
(^ 83 / 2 ) and to Eground (^Ii 5 / 2 ) transitions 

(Fig. 6 a). This means that for a given laser power, rel¬ 
ative intensities will vary from particle to particle (Fig. 
6 b). It is, therefore, improbable that a comparison of 
the individual particles’ emission ratio of the two tran¬ 
sitions to a calibrated ensemble would yield an accurate 
temperature. For this reason, we use the Boltzmann dis¬ 
tribution of these electronic states solely in a qualitative 
manner to confirm the trend of a decreasing or increasing 
temperature. 
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FIG. 7. DDSCAT Simulations of Optically Trapped 

YLF. a, Optically-trapped YLF particle in an aqueous fluid 
chamber. Angles are used to describe numerical discrete 
dipole approximation calculations in panels c-f. b, Cross sec¬ 
tional plot of the internal electric field magnitude within a 
YLF particle with dimensions 7.6 fim & 3.9 /rm at a trap¬ 
ping wavelength of 1020nm. c, Radiation pressure on a YLF 
particle shown in b. d, Torque calculations as a function of 
angular coordinates for the YLF particle shown in b. e, Radi¬ 
ation pressure on a 200 nm long YLF. f, Torque calculations 
as a function of angular coordinates for a 200 nm long YLF 
particle. 


D. Appendix IV: Yb®+ Ion Spacing 

According to the XRD data presented in Fig. Ic and 
detailed in Materials & Methods, the lattice parameters 


of 10%Yb3+:YLF are a = 5.1641 A and c = 10.7177 A 
and the volume of the corresponding unit cell is Vq = 
a^*c = 2.85*10“^® m®. For 1 m® YLF, the number of 
unit cells is n = 1 m®/Vo = 3.5*10^^ m“®. For each unit 
cell, there are 4 rare earth atoms. The total number of 
rare earth atoms in 1 m® is n^ = 4n = 1.4*10^® m“®. For 
10%Yb®“'':YLF, the number of Yb atoms in 1 m® is nyt = 
0.1*n( = 1.4*10^^ m“®. The linear density of Yb atoms 
is p = (nyf,)*(l/3) = 1.11*10® m“^ and the average ion 
spacing is approximately 1 = 1/p = 8.9 A. 

E. Appendix V: Electromagnetic Simulations of 
Optically Trapped YLF 


Discrete Dipole Approximation software DDSCAT 
7.3 [31] was used to simulate laser trap interactions with 
ideal (non-truncated) YLF tetragonal bipyramids in wa¬ 
ter. For a bipyramid of length of 7.6 /rm and base of 
3.9 fiui (Fig. 7b), a dipole density of 2267 dipoles//rm® 
(85,901 dipoles/octahedron) was used. At 1020 nm, re¬ 
fractive indices of YLF and water are 1.44831 and 1.327, 
respectively. Although YLF is a birefringent crystal, this 
has been ignored in the simulation based on the small 
birefringence (~0.02), which will only be reduced by the 
lattice mismatch of the many domains within a single 
octahedron. Also, the approximation of a laser trap by 
a plane wave becomes less appropriate as the size of the 
particle relative to the beam spot gets large. By varying 
0, $, and /3 (see Fig. 7a) and calculating the radiation 
pressure and torque we are able to predict the optimal 
orientation for trapped particles since we assume the par¬ 
ticle will orient so as to minimize these forces and maxi¬ 
mize the trapping potential. One intriguing result of the 
simulations is the particle size dependence of orientation. 
For example, a tetragonal bipyramid with a length of 7.6 
/im is expected to align with the long axis perpendicu¬ 
lar to both k and E, or 0=90°, $=90°, and /3=0° (Fig. 
7c/d), whereas a similar particle with a length of 200 nm 
is predicted to align its long axis to be parallel with the 
direction of propagation (Fig. 7e/f). It is also notable 
that the particle is able generate significant internal and 
near fields as demonstrated in Fig. 7b. 


















